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ABSTRACT 


Experimental  studies  demonstrate  a  significant  effect  of  specimen  size,  a  /  W  ratio  and  prior 
ductile  tearing  on  deavage  fiiacture  toughness  values  (jJe)  measured  in  the  ductile-4o- 
brittle  transition  region  of  ferritic  materials.  In  the  lower-transition  region,  deavage  frac¬ 
ture  often  occurs  under  conditions  offarge  scale  yielding  but  without  prior  ductile  crack 
extension.  The  increased  toughness  develops  when  plastic  zones  formed  at  the  crack  tip  in- 
tonct  with  neaxi^  specimen  surfaces  which  relaxes  crack-tip  constraint  (stress  triaxial- 
ity).  In  the  mid-supper  transition  region,  small  amounts  of  ductile  crack  extension 
(often  <1-2  mm)  routinely  precede  termination  of  the  J-Aa  curve  by  brittle  firacture. 
Large-scale  yielding,  coupled  with  small  amounts  of  ductile  tearing,  magnifies  the  impact 
of  small  variations  in  microscale  material  properties  on  the  macroscopic  firactiire  toughness 
which  ctmtributes  to  the  large  amount  scattm*  observed  in  measured  c7c~values. 

Previous  work  by  the  authors  described  a  micromedianics  fracture  model  to  correct  mea¬ 
sured  «4-values  for  the  mechanistic  effects  of  large-scale  yielding.  This  new  work  extends 
the  modd  to  also  indude  the  influence  of  ductile  crack  extension  prior  to  deavage.  Ductile 
crack  extensions  of  10-15  x  the  crack-tip  opening  displacement  at  initiation  are  considered 
in  plane-strain,  finite  element  computations.  The  finite  element  results  demonstrate  a  sig¬ 
nificant  elevation  in  crack-tip  constraint  due  to  macroscopic  “sharpening”  of  the  extending 
tip  relative  to  the  blirnt  tip  at  initiation  of  growth.  However,  this  effect  is  offset  partially 
by  the  additional  plastic  deformation  assodated  with  the  increased  applied  J  required  to 
grow  the  crack.  Ihe  initial  a/W  ratio,  tearing  modulus,  strain  hardening  exponent  and 
specimen  aua  interact  in  a  complex  mann^  to  define  the  evolving  near-tip  conditions  for 
deavage  fracture.  The  paper  explores  development  of  the  new  model,  provides  necessary 
graphs  and  procedures  for  its  application  and  demonstrates  the  effects  of  the  model  on  firac- 
ture  data  sets  for  two  pressure  vessel  steels  (A533B  and  A515). 
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Effects  of  Prior  Ductile  Tearing  on  Cleavage 
Fracture  Toughness  in  the  Transition  Region 

1.  INTRODUCTION 

Experimental  studies  consistently  reveal  large  scatter  in  the  measured  values  of  cleavage 
fracture  toughness  for  ferritic  steels  tested  in  the  ductile-to-brittle  temperature  transition  re¬ 
gion  (see  for  example:  12, 13, 18, 28, 30, 34-37].  Moreover,  data  sets  of  seemingly  identical  spec¬ 
imens  frequently  show  a  mix  of  cleavage  fractures  with  and  without  small  amounts  of  prior 
stable  tearing;  firacture  toughness  values  for  specimens  with  prior  growth  often  are  significant¬ 
ly  larger  than  tou^mess  values  for  specimens  without  growth.  The  large  amount  of  scatter, 
coupled  with  specimens  that  exhibit  tearing  prior  to  cleavage  and  others  that  show  no  tearing, 
greatly  complicates  the  interpretation  of  toughness  data  to  define  meaningful  values  for  ap¬ 
plication  in  fracture  assessments  of  structures. 

Scatter  in  fracture  toughness  values  (Je)  observed  in  the  transition  region  arises  frx>m  the 
complex  interplay  between  the  mechanistic  effects  of  cradc-tip  constraint  (stress  triaxiality) 
and  the  stochastic  variability  of  cleavage  resistance  at  the  microstructural  level.  Interaction 
of  crack-tip  plastic  zones  wi&  overall  plastic  deformation  of  the  test  specimen  creates  a  highly 
nonlinear  relationship  between  J  and  near-tip  stress  triaxiality.  For  specimens  experiencing 
significant  levels  of  plastic  deformation,  large  increases  of  J  must  be  imposed  to  produce  micro- 
structurally  significant  increases  in  the  near-tip  stress  fields.  Plasticity  remote  from  the  crack 
shields  the  near-tip  region  frnm  the  full  effect  of  the  far-field  J,  which  scales  with  the  global 
plastic  deformation.  This  shielding  phenomenon  magnifies  the  impact  of  small  variations  in 
the  microstructural  cleavage  resistance,  and  contributes  a  large  measure  of  the  scatter  ob¬ 
served  in  fracture  data.  Dodds  and  Anderson  [3, 7]  proposed  a  st^sed-volume,  micromechan¬ 
ics  model  for  cleavage  fracture  to  correct  measur^  (^-values  for  the  medianistic  effects  of 
large-scale  yielding.  The  model  defines  a  constraint  "corrected”  Jo-value,  denoted  JssY 
lier  work,  for  each  measured  <^-value  at  fracture.  The  corrected  Jo-value  generates  the  same 
stressed-volume  of  material  at  the  crack  tip  in  small-scale  yielding  (SSY  with  7^0)  that  Jc 
generates  at  the  crack  tip  in  the  test  specimen  at  fracture.  Under  increased  loading  in  test  spec¬ 
imens,  the  ratio  of  applied  J  to  Jq  ( J/ Jq)  exceeds  unity  as  large-scale  yielding  causes  constraint 
to  fall  relative  to  SSY  reference  condition.  For  example,  the  J  /Jo  ratio  grows  much  larger  than 
4-5  in  shallow  notdied  bend  specimens  (a/W<0.2). 

When  small  amounts  of  ductile  tearing  precede  failure  by  cleavage,  this  picture  becomes 
more  complex.  Stress  fields  ahead  of  a  growing  crack  differ  fium  those  of  a  stationary  track  and 
growth  exposes  additional  volmnes  of  material  to  high  stresses  as  the  tip  moves  forward.  Wal¬ 
lin  [35]  proposed  simple  extensions  of  his  WST  statistical  model  [34]  to  include  the  first  order 
effects  of  stress  field  changes  and  additional  volume.  Recent  work  by  Varies  and  Shih  [32]  on 
steady-state  crack  growth  in  SSY  and  by  ODowd,  et  al.  [22]  for  transient  growth  in  SSY  dem¬ 
onstrates  clearly  the  development  of  high  constraint  conditions  ahead  of  the  growing  crack, 
even  when  growth  b^ins  under  low  constraint  conditions  (T/ao<0).  These  finite  element 
analyses  reveal  a  significant  elevation  of  stress  triaxiality  caused  by  macroscopic  sharpening 
of  the  extending  tip  and  by  additional  elastic  strains  as  the  cradc  extends  forward  into  pre¬ 
viously  imyielded  material.  La  frracture  test  specimens,  however,  these  effects  may  be  partially 
or  completdy  ofibet  ly  the  large  amount  of  plastic  deformation  present  at  Jjc  and  the  additional 
plastic  defonnation  associated  with  the  increased  applied  J  required  to  grow  the  crack.  The 
initial  a/W ratio,  the  tearing  modulus,  strain  hardening  exponent  and  specimen  size  interact 
in  a  complex  manner  to  define  the  evolving  near-tip  conditions  for  cleavage  fii^cture. 


1 


This  work  extends  the  stressed-volume  constraint  model  for  cleavage  firacture  to  include 
the  mechanistic  effects  of  small  amounts  of  ductile  tearing  on  the  near-tip  stress  fields.  Plane- 
strain,  finite  element  analyses  are  described  for  SECB)  specimens  having  both  shallow  and 
deep  cradcs  (a  /  W  s  0. 1, 0.5)  for  materials  with  hardening  exponents  n  s  5, 10.  The  models  are 
loaded  to  prescribed  levels  of  J  =  J/c  after  which  crack  growth  occurs  consistent  with  the  speci¬ 
fied  resistance  curve,  -  Aa.  A  parameter  study  for  a  range  of  Jfe  and  tearing  modulus  values 
provides  the  stress  fiel^  ahead  of  the  growing  cracks  needed  to  construct  J/Jo  ratios  of  the 
constraint  correction  model.  In  this  process^  .we  assume  that  ductile  crack  growth  prior  to  cleav¬ 
age  firacture  introduces  additional  complexity  only  in  the  resolution  of  near-tip  stresses,  i.e., 
that  the  microstructural  resistance  to  cleavage  initiation  remains  stress-controlled  and  unaf¬ 
fected  by  the  small  amounts  of  prior  stable  tearing  (usually  1-2  mm  or  less). 

Given  a  measured  touf^ess  value  at  deavage  (Jc)  and  estimated  values  for  the  tearing 
modulus  and  Jjc ,  the  new  model  corrects  for  the  combined  effects  large-scale  yielding  and  duc¬ 
tile  crack  extension.  Necessary  graphs  and  procedures  are  provided  in  the  last  section  to  apply 
the  model  to  experimental  data  sets.  Statistical  analyses,  e.g.,  Weibull  methods,  developed  for 
interpretation  of  deavage  firacture  data  should  be  applied  to  data  sets  of  corrected  Jo-^alues, 
with  no  distinction  then  needed  between  corrected  values  for  specimens  with/without  prior 
tearing.  The  paper  condudes  with  example  applications  of  the  proposed  model  to  data  sets  for 
two  pressure  vessel  steels  (A515,  A533B). 
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2.  MICROMECHANICAL  CONSTRAINT  CORRECTIONS 


Dodds  and  Anderson  [3, 7]  show  that  by  quantifying  the  effects  of  large-scale  yielding  (LSY) 
on  the  relationship  between  microscale  crack  driving  force  (e.g.  near-tip  stresses  and  strains) 
and  macroscale  crack  driving  force  (e.g.  J,  CTOD),  the  in-plane  size  effect  on  cleavage  fracture 
toughness  can  be  predicted  rigorously  without  resort  to  empirical  arguments.  Statistical  anal¬ 
yses  should  be  applied  only  to  constraint  "corrected”  measures  of  fractiire  toughness  [1,37\. 
Size  effects  become  steadily  more  pronoxiqp^  as  load  increases  due  to  the  deviation  of  crack- 
tip  region  deformations  frmn  the  small-sccde  yielding  (SSY)  conditions  iinder  which  single  pa¬ 
rameter  fracture  mechanics  apply  123,27}.  When  J (or  equivalently  CTOD)  no  longer  describes 
the  crack-tip  fields,  a  microme<^anics  failure  criteria  is  required  to  establish  the  near-tip 
conditions  at  fracture.  Finite  element  anal]rsis  provides  a  means  to  quantify  the  geometry  de¬ 
pendent  relations  between  these  micromechanical  failure  conditions  and  macroscale  crack 
driving  force.  This  permits  (in  principle)  prediction  of  fracture  in  any  body  firom  toughness  val¬ 
ues  measured  using  standard  specimens.  This  section  describes  the  micromechanism  for  cleav¬ 
age  fracture  developed  in  previous  work  to  correct  measured  fracture  data  for  constraint  loss 
due  to  large-scale  yielding  and  then  discusses  the  basis  for  extending  the  methodology  to  in¬ 
clude  small  amounts  of  ductile  tearing  prior  to  cleavage  fracture. 

2.1  Transgranular  Cleavage  Mechanism 

A  number  of  micromechanical  models  for  transgranular  deavage  fracture  have  been  pro¬ 
posed,  most  derive  firom  weakest-link  statistics  (see  review  by  Wallin  [36]).  The  weakest-link 
models  assume  the  largest  or  most  favorably  oriented  fracture-triggering  partide  controls  the 
deavage  failure.  The  actual  trigger  event  involves  a  local  GrifBth  instability  of  a  microcrack 
which  forms  at  a  microstructural  feature  such  as  a  carbide  or  indusion;  satisfaction  of  of  the 
Griffith  energy  balance  occurs  when  the  critical  stress  is  reached  in  the  vicinity  of  the 
microcrack  [2^.  The  size  and  location  of  the  triggering  microstructural  feature(s)  dictate  the 
fracture  toughness  and  contribute  to  the  scatter  observed  in  results  of  deavage  fracture  tests. 

The  Griffith  instability  criterion  implies  fracture  at  a  critical  normal  stress  near  the  crack 
tip;  the  statistical  sampling  aspect  of  the  mechanism  (i.e.,  the  probability  of  finding  a  triggering 
microfeature  near  the  crack  tip)  suggests  a  dominant  role  for  the  volume  of  material  within  a 
process-zone  over  which  the  opening  mode  stress  exceeds  a  threshold  value  suffident  to  initi¬ 
ate  cleavage.  The  probability  of  deavage  fracture  in  a  cracked  specimen  may  then  be  expressed 
in  the  following  general  form: 

F  =  FlV(ai)l  (1) 

where  F  is  the  failure  probability.  is  the  maximum  prindpal  stress  at  a  point  and  V(ai)  is 
the  volume  sampled  over  whidi  the  prindpal  stress  is  equal  to  or  greater  than  ai.  This  form 
of  F  applies  to  any  fracture  process  controlled  by  maximum  prindpal  stress,  not  just  weakest- 
link  failure  [2,37].  In  particular,  the  F  criterion  of  Eq  1  does  not  require  material-spedfic  as¬ 
sumptions  for  the  distribution  and  strength  of  cleavage  triggering  particles. 

The  present  methodology  does  not  attempt  to  predict  absolute  values  of  Jc  Grom  metallurgi¬ 
cal  parameters  that  descri^  the  distribution  and  strength  of  cleavage  triggering  particles. 
Rather,  the  micromechanical  model  predicts  the  variation  of  fracture  toughness  with 
constraint  changes  for  a  given  material/temperature  by  scaling  to  a  reference  constraint  condi¬ 
tion.  The  crack-tip  stress  fields  in  a  test  specimen  are  compared  to  those  for  SSY  with  T=0.  A 
J-like  parameter,  denoted  Jq,  is  obtained  firom  this  comparison  to  the  reference  solution.  Jo  is 
the  J  to  which  the  SSY  model  (infinite  body)  must  be  loaded  to  achieve  the  same  stressed  vol¬ 
ume,  and  thereby  the  same  likelihood,  F,  of  cleavage  fracture,  as  m  a  finite  body.  With  this  mod- 
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el,  the  assurance  of  equally  stressed  volumes  of  material  at  firacture  does  not  imply  equivalence 
of  the  individual  stress  components  beyond  SSY  conditions,  however,  such  a  requirement  does 
not  appear  strictly  essential. 

2.2  Constraint  Corrections 

Consider  the  family  of  near-tip  stress  fields  applicable  to  finite  bodies  described  in  terms 
of  the  J-Q  theory  [20,21] 


where  Qoq  denotes  an  approximately  constant,  hydrostatic  stress  in  the  forward  sector, 
1 0 1  <  90°,  that  defines  the  local  deviation  of  the  finite  body  stresses  firom  those  of  the  small- 
scade  yielding  reference  field  (Q  =  0).  By  employing  the  above  family  of  near-tip  stress  fields, 
the  maximum  principal  stress  also  has  the  form 


For  any  given  value  of  Q  and  6,  oiIoq  decreases  monotonically  once  r  extends  beyond  the  fin¬ 
itely  deformed  region  of  r  2  J/oq.  Rearrangement  of  the  above  expression  furnishes  a  relation 
for  the  distance  r  as  a  function  of  6  and  a\laQ  as 

r  =  ^giWia^/aQ,  Q)  .  (4) 

Consider  a  particular  level  of  the  principal  stress  ai/ao-  The  area  A  over  whidr  the  principal 
stress  exceeds  oi/cq  is  given  by 

A  =  Q),  ^  =  5  I Oi/oof  '  ^5) 

The  area  enclosed  by  the  contour  of  level  oi/ao  depends  on  J  as  well  as  the  triaxiality  of  the 
near-tip  fields  identified  with  Q.  lb  fix  ideas,  let  Aq  and  Jq  designate  the  area  and  J  associated 
with  the  Q  =  0  field,  and  let  Aj^  and  Jfb  designate  the  area  and  J  associated  with  a  crack  in 
a  finite  body  with  Q^O.  Tlien  we  have 

X 

Aq  =  ^0*^1  g^&lOi/oQ,  Q  =  0)dd  (6) 

^  ~X 

and 

X 

Ap2  ~  ^FB  ~  ^  j  g^^i 

Upon  initial  loading  of  the  finite  body,  Q  =  0  so  that  hjrs  =  ho-  As  the  load  increases,  plasticity 
spreads  over  the  body,  Q  becomes  non-zero,  and  ApB  begins  to  deviate  firom  ho. 

For  a  given  material  and  temperature,  firacture  occurs  upon  the  attainment  of  equivalent 
stressed  volumes  (Afs  ^  thickness  B)  for  cleavage  firacture  in  different  specimens.  The  ratio 
of  applied  J-values  in  a  finite  body  and  the  reference  Q  =  0  stress  state  that  generate  equiva¬ 
lent  stressed  volumes  is  fo\md  by  equating  areas  in  Eqs  6  and  7  to  yield 
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The  J  ratios  are  evaluated  iising  Eq  8  at  each  loading  level  and  for  a  range  of  principal  stress 
values.  The  ratio  quemtifies  the  size  and  geometry  dependence  of  cleavage  fracture  toughness. 
Consider,  for  example,  a  test  specimen  that  fails  at  Jc  =  200  kJ/m^.  Suppose  the  computed  ratio 
=  2  at  firacture  (Jps  s  Jc)  in  the  test  specimen;  then  a  much  larger  specimen  made  from 
the  same  material  and  tested  at  the  same  temperature  is  predicted  to  fail  at  Jc  =  100  J/m^  (the 
lager  specimen  fails  under  SSY  conditions).  Similarly,  the  fracture  toughness  ratios  for  test 
specimens  with  the  same  absolute  size  but  varying  crack  depths-to-spedmen  widths,  a/W, 
may  be  quantified.  The  model  predicts  a  sjkarp  increase  in  fracture  toughness  with  decreasing 
a/W  ratio. 

An  often  overlooked  benefit  of  both  the  J-Q  description  of  near-tip  fields  and  this 
constraint  model  for  cleavage  fracture  involves  the  capability  to  model  arbitrary  material 
hardening.  A  finite  element  analysis  of  the  SSY  reference  condition  may  be  executed  using  a 
description  of  the  uniaxial  material  response  other  than  Ramberg-Osgood  or  the  linear,  pow¬ 
er-law  models  (finite  strains  and  J2  flow  theory  plasticity  may  be  modeled  as  well).  Fracture 
tert  specimens  (and  structural  components)  are  then  modeled  using  the  same  procedures  with 
Q  and  the  constraint  model  for  cleavage  constructed  as  above.  These  approaches  rely  on  the 
concept  of  similarity  under  SSY  conditions;  a  concept  much  broader  than  is  inferred  by  asymp¬ 
totic  analyses  which  invoke  incompressibility,  power-law  hardening  and  small-strain  theory 
to  develop  expressions  for  crack-tip  fields. 

The  key  question  to  resolve  with  the  constraint  model  concerns  the  sensitivity  of  the  A^/Ao 
and  Jfs/Jb  ratios  on  oi/oq.  Numerical  computations  [3,14,15]  reveal  a  weak  dependence  on 
ai/oQ  in  common  firacture  specimens  until  ^e  overall  bending  deformation  impinges  signifi¬ 
cantly  on  the  cra(k-tip  field.  Ihis  occurs  more  readily  for  low  amoimts  of  strain  hardening  and 
for  deeply  notched  specimens  and  much  less  so  for  tension  dominated  geometries  and  for  mate¬ 
rials  with  greater  amounts  of  strain  hardening.  The  weak  dependence  of  Jfb  U 0  ratios  on  o\  /  ao 
derives  finm  the  observation  that  stress  fields  among  firacture  specimens  vary  primarily  due 
to  the  hydrostatic  contribution  (as  described  by  the  J-Q  theoiy ).  The  differences  in  hydrostatic 
stress  relative  to  the  SSY  (Q=0)  field  alter  pr^ominantly  the  absolute  size  of  principal  stress 
contours  but  not  the  angular  form  defining  their  shapes  [8). 

2.3  Application  of  the  Constraint  Corrections  in  Fracture  Testing 

Figure  1  provides  typical  results  obtained  through  finite  element  modeling  (plane  strain). 
The  specimens  are  SE(B)8  with  alW  s  0.1, 0.5  modeled  with  a  Ramberg-Osgood  material  hav¬ 
ing  strain  hardening  exponents  n  s  5, 10.  Values  of  Jjtb  and  Jo  are  plotted  on  separate  axes  to 
facilitate  the  constraint  correction  of  experimental  data.  Points  on  the  curves  describe  {JpEt  Jo) 
pairs  that  produce  equal  stressed  volumes  of  material  in  the  test  specimen  and  in  SSY.  Path 
A-B-C  indicated  in  fig.  1  (b)  illustrates  the  procedure  to  correct  for  the  effects  of  large-scale 
yielding  in  measured  cleavage  firactiire  toughness  data  (Jc value  at  A)  by  determining  the  ge¬ 
ometry  independent  fi:acture  toughness  (Jq  value  at  C).  Tb  simplify  application  of  the 
constraint  correction  procedures,  curve  fits  to  finite  element  results  of  the  form  shown  in  Fig. 
1  are  available  for  a  wide  range  of  hardening  exponents  and  a/W  ratios  [8].  Moreover,  crack 
opening  stresses  on  the  remaining  ligament  may  be  used  to  compute  the  Jfb/Jq  ratios  rather 
than  the  more  complex  contour  area  procedure  (81.  Applications  of  this  technique  to  structural 
steels,  including  A36,  A515  and  A533B  [13,30,37],  have  been  successfiil  in  removing  the  geom¬ 
etry  dependence  of  Jc  values. 

Ihree-dimensional  analyses  now  underway  [10]  for  the  same  SE(B)  specimens  (with BxB 
cross  sections)  reveal  significant  differences  with  the  plane-strain  correction  model  for  the 
deep-notch  specimen  but  essentially  no  differences  with  plane-strain  results  for  the  shallow- 
not^  specimen.  The  3-D  analyses  show  that  SSY  conditions  for  the  a/W  s  0.5  specimen  are 
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maintained  to  greater  deformation  levels  {aO{JJ  =  50  vs.  150  for  plane-strain)  on  the  mid¬ 
thickness  plane;  at  positions  on  the  crack  front  mid-way  between  the  center-plane  and  outside 
surface,  the  constraint  levels  fall  below  those  predicted  by  the  plane-strain  model.  For  the 
shallow-notch  specimen,  constraint  conditions  on  the  mid-thickness  plane  closely  match  those 
predicted  by  the  plane-strain  model.  Introduction  of  the  finite  thickness,  B,  requires  modifica¬ 
tion  of  the  constraint  correction  model  to  accommodate  the  interaction  of  thickness  with  a/W 
ratio,  W/B  ratio  and  material  strain  hardening.  The  concept  of  an  effective  thickness,  B^ff,  is 
being  explored  to  describe  constraint  across  the  crack  front  relative  to  the  mid-thickness 
constraint  in  a  meaningful  way  for  deavi^e  fracture.  The  relative  effects  of  crack  growth  in 
a  plane-strain  vs.  a  fully  3-D  model  (with  curved  crack  fi:unts  during  growth  allowed)  appear 
to  be  a  very  important,  yet  currently  unresolved  issue  (the  computational  resources  and  model¬ 
ing  complexity  for  growth  analyses  in  3-D  are  much  greater). 

2.4  Extension  to  Include  Effects  of  Ductile  Tearing 

Small  amounts  of  ductile  crack  extension,  typically  less  than  1-2  mm,  preceding  cleavage 
fracture  alter  the  structure  and  severity  of  the  near-tip  stress  fields.  At  the  length-scale  char¬ 
acteristic  of  the  CTOD  at  initiation  of  tearing,  the  extending  tip  develops  a  sharp  opening  pro¬ 
file  which  elevates  the  stress  and  decreases  the  distance  from  the  crack  tip  to  the  peak  stress 
location.  When  SSY  conditions  prevail,  stresses  at  distances  greater  than  J /aofrom  the  tip  in¬ 
crease  as  well  due  to  additional  elastic  strains  as  the  plastic  zone  pushes  forward  into  unyielded 
material.  Crack  growth  also  subjects  a  larger  volume  of  material,  AV  ec(J/ao)xAax  B,  to  high 
stress  levels  near  the  tip  and  th\is  raises  the  probability  of  sampling  a  critical,  fracture-trigger¬ 
ing  partide. 

At  the  metallurgical  length-scale  characteristic  of  carbides  and  indusions,  the  potential 
exists  for  tearing  and  large  plastic  deformation  to  affect  the  fundamental  cohesive  strength 


»^Fb/(^//ow)  '^Fb/ inflow) 

(a)  (b) 


FIG.  1 — Cleavage  cojistraint  models  for  stationary  cracks  in  SE(B)  specimens  with 
a/WsO.!,  0.5;  n^5, 10;  BIoq  =  500;  using  contour  areas  for  aj  Ioq  =  2.7. 
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FIG.  2 — Schematic  illustration  of  growth  effects  on  constraint  correction  model. 


across  cleavage  planes.  More  likely,  tearing  makes  the  microstructural  crack  front  consider¬ 
ably  more  torturoxis  than  the  blunted  tip  of  the  stationaxy  cradi:,  thus  increasing  the  cleavage 
energy  needed  to  initiate  a  macroscopic  fracture.  In  addition,  material  damage  caused  by  the 
formation  and  growth  of  voids  within  the  zone  of  finite  strains  may  alter  somewhat  the  macro¬ 
scopic  stress  field  [16]. 

lb  extend  the  constraint  correction  model,  these  metallurgical  and  damage  aspects  of  tear¬ 
ing  are  considered  to  have  secondary  importance  rdative  to  the  increased  near-tip  stresses 
and  the  additional  volume  of  material  affected  by  crack  growth.  The  effects  of  ductile  tearing 
thus  appear  as  changes  in  the  stressed  volumes  of  material,  Afb  where  crack  growth 
changes  the  area  enclosed  within  principal  stress  contours  relative  to  a  stationary  cradr  loaded 
at  the  same«^. 


Figure  2  illustrates  schematically  the  impact  of  tearing  on  the  constraint  correction  model . 
The  specimen  follows  the  model  for  a  stationary  crack  imtil  the  onset  of  ductile  tearing  at 
JpB  =  «^/c'  Under  continued  loading,  the  new  model  including  growth  predicts  an  increase  of 
c/o-values  rdative  to  continued  loading  of  the  stationary  cradr,  i.e.,  tearing  increases  the 
constraint.  The  increase  of  Jo-values  due  to  tearing  depend  on  tire  amount  of  plastic  deforma¬ 
tion  in  the  specimen  at  Jjg  and  on  the  ff-curve  followed  prior  to  deavage  firacture.  lb  simplify 
initial  development  of  the  model,  a  constant  tearing  modulus,  Tj  «  (Mfa^dJIda,  is  adopted 
to  describe  the  ff-curve  for  typical  extensions  of  1-2  mm.  Consequently,  the  evolution  of  crack- 
tip  corrstraint  in  simple  fracture  specimens  with  ductile  growth  may  be  expressed  by  a  function 
of  the  form 


JjB 

Jo 


^*^0  Ac  \ 
W  Jj^  ’ 


(9) 


where  C  denotes  the  “constraint”  function  for  deavage  fracture  with  prior  ductile  tearing.  The 
non-dimensional  groups  Wbg/Jj^,  Oq/W,  Aa/W  and  Tj  describe  the  dependence  on  crack 
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growth  and  absolute  size,  while  EIoq  and  n  define  the  material  yield  and  hardening  properties. 
Although  included  Eq  9,  a  weak  dependence  of  C  on  ai/oois  expected. 

Jo-values  derived  firom  both  instantaneous  and  cumulative  contour  areas  are  provided  for 
initisd  development  of  the  growth  correction  model  (see  Fig.  3).  The  instantaneous  areas  en¬ 
closed  within  principal  stress  contoiirs  are  computed  firom  a  snapshot  of  the  stress  field  ahead 
ofthe  growing  crack— points  on  such  contours  all  lie  within  the  forward  sector,  |9|  <90°  (see 
area?  in  Fig.  3).  Figure  3  also  illustrates  the  construction  of  cumulative  contour  areas  using 
the  envelope  of  tnaTimnw  principal  stresses  experienced  by  material  along  the  cra^  plane  dur¬ 
ing  growth.  The  cumulative  areas,  defined  as  areas  enclosed  within  the  contour  envelope,  may 
be  substantially  larger  than  instantaneous  areas  which  leads  to  the  larger  Jo-values  indicated 
in  Fig.  2  for  the  cumulative  area  correction. 

The  proposed  constraint  model  induding  crack  growth  maintains  the  relative  simplidty  of 
the  original  version  whidi  corrects  only  for  (inr-plane)  efiTects  of  large-scale  yielding.  Further¬ 
more,  development  of  the  new  model  to  account  only  for  tearing  effects  on  the  microscale  driv¬ 
ing  force  (the  local  crack  tip  stresses  and  stressed  volume),  implies  that  corrected  toughness 
values,  with  and  without  ductile  tearing  prior  to  deavage,  should  receive  equal  weight  in  sta¬ 
tistical  analyses,  e.g.,  Weibull  models. 


Contour  Envelope  Area  for  ai/ao=A 
OndudesA&^ 


Contour  Area  (A)  for 
oi/oo=A  @  Aa=0 


Instantaneous  Contour 
Area  (A)  for  o-^Ioq-X  @  Aa  >  0 


FIG.  3 — Schematic  illustrating  construction  of  principal  stress  envelope  for  growing 
cracks. 
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3.  NUMERICAL  MODELS  FOR  GROWTH  ANALYSES 


Finite  element  models  of  the  type  shown  in  Fig.  4  are  employed  to  perform  plane  strain, 
nonlinear  analyses  of  SE(B)  specimens  including  small  amoimts  of  ductile  crack  growth.  Both 
shallow  crack  (o/WsO.l)  and  deep  crack  (a/WaO.5)  specimens  are  analyzed.  Symmetry  about 
the  crack  plane  permits  modeling  of  only  one-half  of  each  SE(B)  specimen.  The  finite  elements 
are  conventional  eight-noded  isoparametrics  with  reduced  (2x2)  Gauss  integration.  Square 
elements  in  the  crack-tip  region  and  aloi^  the  crack  plane  are  defined  to  simplify  the  growth 
analyses  and  to  permit  uniform  incremmts  of  crack  extension.  Crack  opening  p'  es  and 
stresses  ahead  of  the  growing  crack  seem  to  be  resolved  adequately  when  the  crack-  .ment 

size, is  smaller  than  •■0.5  x  theCTODatinitiationCd/e).  Stationary  crack  solution,  or  com¬ 
parison  are  also  generated  with  these  models. 

The  finite  element  models  are  loaded  by  displacement  increments  imposed  on  a  single  cen- 
terplane  node  as  shown  in  Fig.  4.  An  improved  form  of  the  traditional  node  release  technique 
is  employed  to  implement  crack  extension.  The  imposed  load-point  displacement  is  adjusted 
iteratively,  concurrent  with  a  one-element  crack  extension,  sudr  that  the  crack  growth  criten- 
on  is  always  satisfied  at  completion  of  the  node  release.  Hus  procedure  retains  the  computa- 


,  X  Rnite  element  mesh  SE(B)  specimen  with  deep  crack 
W  (a/W  =  0.5).  Elements  are  8-node  isoparametrics. 

Le/W  =  0.001 


Crack-plane  elements  are  square  to  simplify  crack 
(b)  growth  analyses.  Lg  denotes  side  length  of  the  square 
crack  plane  elements. 

FIG.  4 — Finite  element  models  for  SE(B)  specimens. 
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tional  convenience  of  the  node  release  technique  while  eliminating  convergence  problems  due 
to  artificial  unloadings  encountered  with  the  traditional  implementation  (see  [d]  for  additional 
details). 

The  uniaxial  stress-strain  rdationship  follows  the  Ramberg-O^ood  model,  e/cg 
a/ug  +  a(a/ag)''.TWotypicalvaluesofthematmialhardeningcoefficient(/i  s  Sand  10)are  se¬ 
lect^  in  the  analyses  of  SE(B)  specimens  to  represent  hi^  (n  »  5)  and  medium(n  s  10)  hard¬ 
ening  materials.  All  computations  use  EIoq  s  500,  Poisson’s  ratio  v  s  0.3  and  a  s  l. 

The  computations  are  carried  out  using  a  large-rotation,  finite  strain  plasticity  model 
derived  firom  incremental  (J2)  theory  A  companion  study  using  a  conventional,  small- 
strain  plasticity  formulation  predicts  essentiaUy  identical  results  for  crack  opening  profiles 
and  stresses  ahead  of  the  growing  tip.  Recent  wcvk  by  Liu  and  Drugan  reach  similar  conclu¬ 
sions  [24\. 

3.1  Crack  Growth  Parameters 

The  specified  value  of  J/e,  relative  to  specimen  geometry  and  material  flow  properties,  sets 
the  severity  of  plastic  defmmation  in  an  SE(6)  specimen  at  the  initiation  of  ductile  growth.  Val¬ 
ues  ciWao/Jje  -  480-»80  are  defined  to  adiieve  conditions  at  initiation  of  growth  ranging  firom 
essentially  SSY  to  fully  formed  plastic  hinges. 

Once  J  reaches  the  specified  initiation  toneless,  Jje,  subsequent  crack  extension  follows 
an  R-curve  defined  by  a  constant  tearing  modulus,  Tj.  Selected  values  for  the  tearing  modulus 
of  15, 40  and  83  represent  relatively  low  valu^  diaracteristic  of  higher^trength  steels  with 
moderate-to-low  strain  hardening.  Values  on  the  order  of  200-400  occur  during  the  early 
stages  of  crack  growth  in  some  lower-strength  steels  with  moderate-to-high  strain  harden¬ 
ing.  These  large  values  approach  one-half  the  slope  of  the  blunting  line,  dJ/ da  s  2ao  (defines 
a  tearing  modulus  of  Tjbl  -  2E/ao).  The  lower  values  of  Tj  specified  for  the  SE(B)  analjrses 
accentuate  the  effects  of  crack  growth  —  mudi  larger  values  of  tearing  modulus  drive  the  re¬ 
sponse  toward  the  stationary  crack  solution,  i.e.,  the  loss  of  constraint  due  to  large  increases 
in  J  and  plastic  deformation  overwhelm  the  increased  constraint  firom  small  cradc  extensions. 
Analyses  for  materials  with  high  tearing  modulus  (200-500)  are  also  performed  to  support  ap¬ 
plication  of  the  growth  constraint  model  to  currently  available  experimental  results. 

3.2  cf-Values  During  Crack  Growth 

J-values  dttring'cradL  growth  analyses  in  the  SE(B)  specimens  are  computed  with  the  do- 
main-integral  scheme  including  only  terms  for  a  nonlinear  elastic  material  model  [29].  Figure 
5  shows  the  typical  variation  of  J-vidues  with  increasing  generalized  radius  of  the  integration 
domain  for  both  a  stationary  crack  and  a  growing  crack.  A  domain  *ring”  denotes  a  sequence 
of  coimected  elements,  enclosing  the  tip  at  a  remote  distance,  whidi  make  non-zero  contribu¬ 
tions  to  J.  Starting  with  the  ring  of  elements  at  the  current  crack  tip,  J  is  computed  indepen¬ 
dently  over  each  successive  ring  of  elements.  Figure  5  (a)  shows  a  weak  path  dependence  of  J 

on  the  dennain  size  (i.e.,  distance  ofringfinm  tip)  for  the  stationary  cra^  loaded  to  'nre 

small  degree  of  path  dependence  just  outside  the  crack-tip  elements  arises  firom  nonpropor¬ 
tional  straining  in  the  incremental  plasticity  theory.  Once  the  domains  extend  beyond  10%  of 
the  remaining  ligament,  Jbecomes  path-independent;  this  value  is  referred  to  here  as  the  far- 
field  J  and  is  denoted  Jf^. 

In  contrast  to  the  stationary  crack,  J-values  for  the  growing  crack  exhibit  a  much  strong&r 
path-dependence  as  shown  in  Fig.  5  (b).  Essentially  zero  J-values  are  computed  over  domams 
very  dose  to  the  cradc-tip.  N^ative  J-values  are  obtained  for  domains  within  about  5%  of  the 
remaining  ligament  (in  this  example)  which  pass  throu^  a  large  portion  of  the  elastic  unload- 
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Domain  Radius  /  (W-a)  Domain  Radius  /  (W-ao) 


(a)  Stationary  Crack  (b)  Growing  Crack 

FIG.  5 — J-integral  variation  over  domain  rings  for  both  stationary  and  growing  cracks 
for  SE(B)  with  a/WsO.5. 


ing  zone  zone  behind  the  growing  tip.  Figure  5  (b)  shows  that  the  J-value  at  Aa  s  10  x  d/e 
readies  a  nearly  path-independent  value  once  domains  extend  b^ond  20-25%  of  the  remain* 
ing  ligament.  The  necessity  for  J  evaluation  at  large  distances  the  growing  crack  limits 

the  maximum  amount  of  crack  growth  permitted  in  the  SE(B)  analyse.  As  die  largest  domain 
is  redefined  for  eadi  increment  of  growth,  the  domain  may  extend  into  the  more  coarse  mesh 
remote  from  the  tip  or  readi  the  specimen  edges  ahead  of  or  behind  the  crack  tip. 

In  the  experimental  determination  of  a  Jg-ha  response,  tf-values  are  estimated  as  a  firac- 
tion  of  elastic  and  plastic  work  done  on  the  specimen.  By  treating  the  finite  dement  P-Au  n 
results  (induding  crack  growth)  as  experimental  records,  a  comparison  of  the  domain  integral 
and  experimental  J  definitions  becomes  possible.  For  deep  notch  SE(B)  specimens,  ASTM 
E1152-87  [4\  outlines  the  procedure  to  estimate  J-values  (Jastm)  experimentally  mea¬ 

sured  load,  load-line  displacement  and  cradc  length  records.  Clhe  ASTM  procedure  assumes 
nonlinear  elastic  material  response  as  well  and  thus  the  Jastm  values  are  often  referred  to  as 
“deformation”  J-values,  Jo).  For  both  hardening  exponents  (n  s  5  and  10),  the  finite  element 
J^  and  Jastm  values  are  indistinguishable  over  the  full  cra^  growth  history.  Sudi  dose 
agreement  of  finite  element  and  experimental  estimates  for  J-values  with  cradc  growth  is  es¬ 
sential  for  subsequent  application  of  the  growth  constraint  modd. 
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4.  GROWTH  EFFECTS  ON  CRACK-TIP  FIELDS 


4.1  Crack  Opening  Profiles 

Figure  6  compares  crack  opening  profiles  for  stationaiy  crack  and  growing  cracks  in  SE(B) 
specimens.  Hie  effects  of  different  tearing  modulus  are  shown  for  crack  growth  initiated  at  the 
same  nonnalized  value  oUjc.  Profiles  for  growing  cracks  exhibit  a  macroscopically  ‘sharp*  tip 
for  all  values  of  tearing  modulus,  in  contrast  to  the  severely  blunted  tip  of  the  stationary  crack 
While  the  different  values  of  tearing  modulus  stron^y  affect  the  immediate  crack  opening 
angle  (CTOA),  as  indicated  on  the  figure,  there  is  negligible  effect  upon  further  opening  bdiind 
the  initial  (blunted)  tip  when  displacements  are  normalized  hyJ/oQ.  Fdlowing  a  short  tran¬ 
sient  period  when  the  crack  grows  through  the  first  element  at  &e  stationaiy  tip,  the  cradc-tip 
opening  angles  remain  constant  and  increase  with  the  specified  Tj.  The  constancy  of  CTOA 
with  Tj  observed  previously  in  SSY  analyses  of  crack  growth  [2d,29]  carries  forward  to  the 
SE(B)  specimens. 

4.2  Crack  Opening  Stresses  nw  T.igi>m«tnt 

Figure  7  illustrates  the  typical  effects  of  growth  on  the  crack  opening  stress  whidi  acts  on 
the  remaining  ligament.  The  analyses  consider  two  geometries  (a/ W  =  0.1, 0.5)  and  two  hard¬ 
ening  exponents  (n  s  5, 10),  with  results  for  each  combination  of  these  parameters  shown  for 
one  value  of  the  tearing  modulus  (7^  =  83).  Crack  growth  begins  in  the  analyses  at  the  common 
value  Jje  =  W(7o/240.  The  results  for  eadi  case  contain  two  curves;  one  for  the  growing  cradi 
of  length  oq  +  Aa  and  one  for  the  stationary  cradc  of  length  oo,  with  both  configurations  loaded 
to  the  same  J.  The  amount  of  crack  extension  relative  to  the  CTOD  at  initiation  varies  with 
the  strain  hardening  even  though  the  final  J-values  are  in  eadi  case.  For  reference,  each 

figure  also  includes  the  corresponding  steady—  growth  solution  in  SSY  given  by  Vaxias  and  Shih 
[32]. 

Sharpening  of  the  crack  tip  durii^  growth  shifts  the  location  of  the  peak  stress  nearer  to 
the  tip,  with  the  largest  (relative)  shift  for  the  lower  amount  of  strain  hardening.  The  peak 


FIG.  6 — Opening  profiles  for  stationary  and  growing  cracks  for  shallow  and  deep  notch 
SE(B)  specimens. 
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stress  location  for  both  high  and  moderate  hardening  converges  to  r«0.25(//ao-  As  expected, 
the  shallow  crack  specimens  reveal  significant  loss  of  constraint  relative  to  the  deep  cradc  spec¬ 
imens,  with  stresses  for  n  s  10  showing  the  most  differences  (compare  stresses  at  the  same 
rl(J /ao)]-  At  distances  firom  the  tip  greater  than  about  J /oq,  crack  growth  in  the  SE(B)  speci¬ 
mens  does  not  alter  significantly  stresses  on  the  remaining  ligament.  This  behavior  contrasts 
sharply  with  stresses  for  the  steadily-growing  crack  in  SSY;  stresses  larger  than  those  for  the 
non-growing  crack  are  predicted  over  much  greater  distances  firom  the  tip  [32].  IntheseSECB) 
specimens,  crack  growth  initiates  and  gnn^s  under  conditions  of  large-«cale  yielding  on  the 
ligament  which  appears  to  suppress  other  mechanisms  of  stress  elevation  beyond  the  sharpen¬ 
ing  effect. 

For  the  n  =  10  material,  the  amplitude  of  peak  stress  during  crack  growth  exceeds  the  am¬ 
plitude  of  the  corresponding  stationary  crack  stress  and,  for  the  deep  notch  specimen,  the  peak 
stress  reaches  the  steady  growth  value  in  SSY.  For  the  Ugh  hardening  material  (ns5),  the 
growth  stresses  appear  to  fall  bdow  the  peak  values  without  growth.  The  relative  coarseness 
of  the  element  mesh  is  clearly  evident  in  the  stress  fields  without  growth  for  the  n  =  5  material; 
the  cra(^-tip  opening  displacement  for  this  case  is  comparable  to  the  tip  element  size.  We  ex¬ 
pect  that  re-analjrsis  using  a  finer  mesh  at  the  tip  forn  =  5  would  predict  somewhat  lower  peak 
stresses  for  the  stationary  crack.  Comparisons  of  crack  opening  stresses  computed  for  other 
cases  reveal  that  variations  of  Jje  and  the  strain  hardening  exponent  affect  the  stresses  more 
significantly  than  variations  of  the  tearing  modulus  when  distances  are  scaled  by  J/oq  as  in  Fig. 
7. 

4.3  Principal  Stress  Contours 

Crack  opening  stresses  on  the  remaining  ligament  may  not  reflect  changes  in  stresses  away 
firom  the  crack  plane  due  to  growth.  Figure  8  (a-b)  show  the  effects  of  crack  growth  on  the 
instantaneous  principal  stress  contours  for  shallow  and  deep  crack  SE(B)  specimens.  Cracdc 
growth  begins  at  the  ccnnmon  value  Jjc  =  Wao/240  with  results  given  in  the  figure  for  7^  =  83 
and  n  =  10.  Comparisons  for  contours  other  tlmnoo/oo  =  2.7  shown  in  the  figure  reveal  similar 
trends.  In  each  case  the  stationary  crack  and  growing  crack  contours  correspond  to  the  same 
applied  J. 

By  using  the  current  value  of  J/oq  to  normalize  distances  firom  the  cradc  tip  on  plots  of  prin¬ 
cipal  stress  contours,  deviations  in  the  stress  triaxiality  due  to  crack  grow^  and  increased 
loading  Garger  J  value)  become  very  dear.  Under  SSY  conditions  the  contoiue  for  a  specified 
ai/ao  ratio  coindde  at  all  tf-values  when  expressed  on  these  normalized  axes.  When  contours 
in  a  specimen  at  different  loads  coindde,  the  near-tip  constraint  remains  xmchanged  (the  abso¬ 
lute  size  of  the  contour  increases  with  <7  but  at  the  same  rate  for  the  two  loading  levds).  When 
constraint  increases  (decreases)  under  increased  loading  or  crack  growth,  the  normalized  con¬ 
tours  must  also  increase  (decrease)  in  size. 

The  normalized  contours  for  the  growing  crack  are  distinctly  larger  than  the  contours  for 
the  stationary  crad^espedally  for  the  shallow  crack  SE(B)  specimen.  For  the  a/ W  =  O.Ssped- 
men,  the  opening  mode  stresses  on  the  remaining  ligament  do  not  indicate  correctly  the  effects 
of  crack  growth  on  stresses  elsewhere  in  the  forward  sector.  The  contours  shown  in  the  figure 
are  those  at  the  maximum  crack  growth  considered  of  10  x 6^;  contours  examined  at  less 
crack  growth  show  somewhat  larger  differences. 

Figure  8  (c-d)  summarize  the  areas  endosed  by  the  instantaneous  contours  for  a  range  of 
ai/ao  ratios.  Areas  are  normalized  by  the  similarity  factor  (in  SSY  with  Aa  s  0,  the  quan¬ 
tity  remains  constant  under  increased  loading).  Three  sets  of  results  for  com¬ 

parison  are  shown  on  each  figure:  (1)  areas  for  SSY  without  crack  growth,  (2)  areas  for  the 
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FIG.  7 — Effects  of  growth  on  crack  opening  stress  in  SE(B)  specimens.  For  each  a/W, 
growing  cracks  and  stationary  cracks  are  loaded  to  the  same  J. 
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FIG.  8 — Effects  of  crack  growth  on  the  instantaneous  principal  stress  contours  in  SE(B) 
specimens.  For  each  a/W,  growing  cracks  and  stationary  cracks  are  loaded  to  the  same  J. 
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SE(B)  spedxnen  at  the  maziinum  amount  of  crack  extension  considered,  and  (3)  areas  for  the 
SE(B)  specimen  without  crack  growth  but  loaded  to  the  same  J-value  as  the  growing  crack. 
Figure  8  (c)  for  the  shallow  crack  demonstrates  the  dramatic  loss  of  constraint  caused  by  the 
extensive  plastic  drformation;  the  areas  differ  by  a  factor  of  100  &om  SSY  for  a\la\  =  2.7. 
Crack  growth  clearly  increases  the  constraint,  elevating  the  areas  by  factors  of  2-3  above  sta¬ 
tionary  crack  levels.  Tlie  deep  crack  specimens  maintain  much  higher  constraint  as  evidenced 
by  the  smaller  differences  in  contour  areas  compared  to  SSY.  Crack  growth  elevates  neaz^tip 
stress  triaxiality  and  increases  areas  enclosed  within  the  principal  stress  contours.  However, 
the  effects  di  growth  on  the  contours  are  much  less  pronounced  than  for  the  shallow  crack. 

Not  surprisingly,  the  differences  in  contour  areas  between  the  stationary  and  growing  crack 
increases  with  the  oi/oq  ratio,  that  is,  contours  near  the  growing  tip  sense  the  increased  peak 
stress  caused  by  the  sharpening.  Hie  trends  in  contour  areas  shown  in  these  figures  are  entire¬ 
ly  consistent  with  observations  about  growth  effects  on  constraint  made  previously  on  the  basis 
of  crack  opening  stresses  acting  on  the  remaining  ligament. 

In  summary,  our  extensive  studies  [31]  of  SE(B)  specimens  with  crack  growth  to  examine 
the  effects  of  h,  Tj  reveal  a  complex  interaction  of  constraint  trade-ofb.  Increased  strain 
hardening  increases  the  level  of  constraint  and  reduces  the  effect  of  growth  on  the  near  tip 
fields.  Larger  values  of J/e  increase  the  impact  of  crack  growth  on  the  near-tip  stresses — crack 
growth  has  more  effect  when  the  stress  fields  at  initiation  have  lost  significant  constraint  rela¬ 
tive  to  SSY.  Larger  values  of  the  tearing  modulus  reduce  the  net  change  in  the  constraint  due 
to  growth — a  large  tearing  modulus  requires  large  J  increases  to  extend  the  crack  whidi  in- 
tum  accelerate  the  constraint  loss  by  additional  large-scale  yielding.  Consequently,  crack 
growth  has  more  effect  as  Tj-»0  and  n^ligible  effect  as  Tj-*Tjbl  (tl^e  blunting  line  slope). 


6.  TEARING  EFFECTS  ON  CONSTRAINT  MODEL 

The  extensive  finite  element  analyses  of  SE(B)  specimens  that  include  the  effects  of  crack 
growth  provide  a  basis  to  extend  the  previously  developed  constraint  model  for  cleavage  firac- 
ture.  The  SSY  stress  fields  (without  crack  growth)  again  define  a  convenient  reference  condi¬ 
tion  to  assess  constraint  effects.  The  Jfb  ratios  are  computed  using  Eq  8  for  two  afW  ratios 
(0.1, 0.5),  for  two  strain  hardening  exponents  (n=5, 10)  and  for  a  wide  range  of  values  for  J/c 
and  tearing  modulus.  The  computations  j>rovide  ratios  based  on  areas  within  the 

instantaneous  principal  stress  contours  ahd  on  cumulative  areas  within  the  envelope  of  con¬ 
tours  constructed  over  the  histozy  of  crack  growth. 

Figures  9  and  10  show  the  typical  variation  of  Jjtb/Jo  ratios  with  principal  stress  contoxir 
found  in  the  parameter  study.  These  figures  illustrate  the  effects  of  a/W  ratio  and  hardening 
exponent  for  a  common  value  of  Jjc  and  Tj,  and  also  compare  the  Jfb/Jq  ratios  for  instanta¬ 
neous  and  cumulative  areas.  The  solid  lines  denote  ratios  at  the  onset  of  crack  growth, 

Jfb  =  «f/c:  eadi  dotted  line  indicates  the  ratios  at  a  fixed  amount  of  crack  growth.  When  the 
JpB^Jo  ratio  increases,  constraint  (stress  triaxiality)  in  the  specimen  decreases. 

Consider  first  the  shallow  crack  results  shown  in  Fig.  9.  The  significant  features  include: 
(1)  high  strain  hardening  (n  s  5)  limits  constraint  loss;  (2)  JpBf^o  ratios  for  crack  growth  are 
less  dependent  on  the  principal  stress  ratio  than  are  the  ratios  for  stationary  cracks  (for  n=10, 
the  stationary  crack  has  lost  similarity  of  the  principal  stress  contours  at  Jjc );  (3)  Jfb/ Jo  ratios 
based  on  cumulative  areas  are  smaller  than  ratios  based  on  instantaneous  areas,  and  (4)  most 
importantly,  ratios  decrease  relative  to  the  stationary  crack  values  with  increased  crack 

growth  and  increased  loading.  This  last  feature  contrasts  sharply  with  the  behavior  of  station¬ 
ary  cracks  which  show  continually  increaring  Jfb/^O  ratios  under  additional  loading. 
Constraint  in  the  shallow  crack  specimen  with  n  =  10  increases  continuously  with  crack 
growth.  The  response  in  Fig.  9(a)  forn  =  5  computed  using  instantaneous  areas  provides  a  par¬ 
ticularly  interesting  result  in  that  the  constraint  remains  essentially  constant — the  increased 
constraint  caused  by  crack  growth  cancels  almost  exactly  the  constraint  loss  due  to  additional 
plastic  deformation  under  continued  loading. 

Consider  next  the  deep  crack  results  shown  in  Fig.  10.  The  overall  trends  remain  similar 
except:  (1)  Jfb/ Jo  ratios  at  the  onset  of  tearing  are  weakly  dependent  on  the  principal  stress 
ratio  for  both  low  and  high  amounts  of  hardening  and  (2)  Jfb/^0  ratios  for  growth  based  on  cu- 
mvilative  area  for  n  s  10  show  a  more  noticeable  dependence  on  the  principal  stress  contour. 
During  crack  growth  with  the  corresponding  increase  in  J,  the  specimen  deformation  increases 
substantially  and  the  bending  field  impinges  on  the  crack-tip  region. 

Figure  11  recasts  these  Jfb/ Jo  ratios  for  oi/oq  -  2.7  into  the  usual  format  of  the  constraint 
model.  Results  of  similar  computations  for  additional  values  of  Jjc  are  given  as  well.  Solid  lines 
denote  the  stationary  crack  response  while  dashed  lines  denote  the  growth  response  in  terms 
of  the  instantaneous  and  cumulative  contour  areas.  Each  curve  for  the  specific  a/W  ratio  termi¬ 
nates  at  the  same  amount  of  crack  gro>vth  in  the  specimen  (Aomox  is  indicated  for  each  set  of 
results).  A  reference  line  indicating  AJq  ==  ^Jfb  is  provided  on  the  figures  to  aid  in  assessing 
the  relative  constraint  changes. 

On  these  and  subsequent  figures,  we  normalize  J-quantities  using  the  flow  stress  ipfutw) 
rather  thim  oq  to  make  the  results  less  sensitive  to  the  actual  estimate  of  the  hardening  expo¬ 
nent.  For  a  Ramberg-Osgood  material  representation  used  in  the  finite  element  computations, 
an  estimate  for  the  flow  stress  is  given  by 
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FIG.  9 — Effects  of  principal  stress  ratio  and  instantaneous  vs.  cumulative  areas  on  Jq  for 
a  sfialloiv-notched  SE(B)  specimen  (alW^O.l). 
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FIG.  10 — Effects  of  principal  stress  ratio  and  instantaneous  vs.  cumulative  areas  on  Jq  for 
a  deep-ncfUhed  SE(B)  specimen  (a/W=0.5). 
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where  N  =  1/n.  The  remaining  ligament  prior  to  crack  growth,  b  =  W-ao,  provides  a  suitable 
specimen  dimension  for  normalizing  J  in  the  constraint  model. 

This  figure  shows  clearly  the  effects  of  crack  growth  on  the  effective  J  for  cleavage  firacture 
(Jq)  for  increasing  Jje  at  a  common  value  of  the  tearing  modulus  (Tj=40).  All  dashed  lines  in 
Fig.  11  lie  above  the  solid  lines,  which  indicate  constraint  increases  with  crack  growth.  Signifi¬ 
cant  observations  firom  this  figure  include:  (1)  the  constraint  model  using  the  cumulative  areas 
always  predicts  larger  values  of  effective  fracture  driving  force  (Jq);  (2)  the  greater  impact  of 
growth  for  the  shallow  crack  specimens  apparent  by  the  mudti  sharper  deviation  of  the 
growth  curves  fiK)m  the  stationary  crack  curve;  (3)  for  this  combination  of  material  parameters, 
the  cumulative  growth  curves  have  slopes  of  1,  i.e.,  increments  of  JpB  applied  to  the  specimen 
after  the  onset  of  growth  are  transmitted  fully  into  increments  of  increased  driving  force  for 
cleavage. 

5.1  Normalizingt/ie  Effects 

Figure  11,  and  others  constructed  for  a  range  of  tearing  modulus  values,  provides  the  mo¬ 
tivation  for  development  of  a  normalization  procedure  to  minimize  the  dependence  of  the  model 
on  values  of  J/^.  The  family  of  growth  constraint  curves  for  different  values  of  the  tearing 
modulus  (at  a  common  Jje)  simply  translates  with  increased  values  for  Jje  while  retaining  a 
nearly  constant  orientation  relative  to  the  gradient  of  the  constraint  curve  for  the  stationary 
crack.  Figure  12  (a-b)  demonstrate  the  effects  of  J/c  on  the  growth  portions  of  the  constraint 
model  when  viewed  in  this  firamework.  In  Fig.  12  (a),  m  represents  the  gradient  of  the  station¬ 
ary-crack  constraint  model  at  Jie,m  =  dJo/dJpB;  ^  denotes  the  arigle  diange  between  the  tan¬ 
gent  line  to  the  stationary  curve  (which  has  slope  m)  and  a  similar  line  drawn  tangent  to  the 
growth  curve  at  the  same  J/e  value.  When  the  J/e  value  increases  from  to  m  decreases 


^FBl{bOfiow)  ^FB^ibOnow) 

(a)  a/IV  =  0.1  (b)  a/IV  =  0.5 

FIG.  11 — Comparison  of  constraint  corrections  for  growth  hosed  on  areas  enclosed  within 
instantaneous  and  cumulative  principal  stress  contours  for  SE(B)  specimens  (Tj-40, 
oi!  00  =  2.7). 
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FIG.  12 — Devdopment  of  approximate  normalization  scheme  to  minimize  dependence  of 
growth  constraint  rnodel  on  absolute  values  ofJjc- 


significantly  from  mi  to  m2  as  indicated  on  the  figure.  Figure  13  provides  values  of  m  for  SE(B) 
specimens  containing  shallow  and  deep  cracks  for  hardening  exponents  n  s  5, 10. 


Figure  12  (b)  shows  growth  curves  defined  at  increasing  values  of  J/e  after  the  applicable 
Jjc  value  is  subtracted  fimm  the  total  values  of  Jfb  and  Jq.  The  shifted  curves  now  show  only 
the  increments  of  J  that  occur  during  growth,  denoted  AJpB  and  Ac7o>  and  defined  b3r; 


AJo  (“) 

where,  for  the  remainder  of  this  section,  all  i/-quantities  are  implicitly  normalized  by  the  appli¬ 
cable  bc^iow\  ropi'asents  the  value  of  Jo  at  Jjrs  ~  J/e*  Using  these  new  quantities,  the  nearly 
equal  values  of  a  at  the  different  values  of  J/e  become  more  apparent.  Tlie  a  values  esdiibit  a 
weak  dependence  on  J/e  for  the  a  /  Ws  0.5  specimen,  but  are  strongly  independent  of  J/e  for  the 
a  /  W  s  0. 1  specimen.  Since  a  remains  relatively  constant,  the  dependence  of  the  growth  curves 
with  J/e  may  be  approximated  as  a  linear  function  of  the  gradients  (m)  for  the  stationary  crack. 


A  simple  procedure  now  defines  the  relationship  between  Jo  v^ues  for  two  different  J/c  val¬ 
ues.  Consider  points  A  and  B  in  Fig.  12  (a)  which  represent  specimens  with  the  same  Tj  and 
same  (normalized)  amo\mt  of  crack  growth  but  with  different  (normalized)  J/e  values  —  con¬ 
figurations  (1)  and  (2)  thus  may  represent  specimens  differing  in  absolute  size  but  having  the 
same  absolute  value  of  J/e-  The  constraint  corrected  values  for  these  two  points  are  labelled 
J^  and  Jf  in  Fig.  12  (a).  During  crack  growth  along  (2)  to  reach  point  B,  JpB  changes  by 
A  J|g  =  A  JjE^;  for  growth  along  (1)  to  reach  A  the  change  in  AJ^  is  also  AJ^  {Tj  and  Aa  are 
the  same).  The  corresponding  changes  in  Jo  values  during  growth  are  shown  as  A  Jg  and  AJ§. 
Using  the  approximation  that  a/  ••02,  Fig.  12  (b)  illustrates  that  the  difference  between  A  J^ 
and  A  Jf  may  be  given  by  the  expression: 

AJ^  -  AJf  -  AJfjf  ’  (mj  -  mj)  . 


m  =  dJoIdJpB  m  =  dJo/dJpB 


dpB/ibotiow)  dpeUbofiow) 

(a)  a/lV  =  0.1  (b)  a/lV  =  0.5 


FIG.  13 — Gradients  (m)  for  stationary  constraint  models  ofSE(B)  specimens  (atW  =  0.1, 
0.5^  n  =  5, 10/  E / ^2  —  SOO). 
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FIG.  14-— GrowrA  portion  of  constraint  model  normalized  to  referertce  point  on  stationary 
crack  solution  for  SE^)  (alW^.l;  n^S,  10). 


Using  this  relationship,  the  constraint  corrected  value  for  point  >1  may  be  expressed  in  terms  of  the 
constraint  corrected  value  for  point  B  as  (solve  above  expression  for  A  Jg) 

jA  ^  jMl)  +  ^jA  .  j/e(l)  +  .  (14) 

Again,all  /  quantities  in  this  expression  are  understood  to  be  nomudized  by  the  applicable  value  of 

For  interpretation  of  experimental  data,  the  corrected  tou^mess  value,  Jq,  is  sought  given 
measured  values  for  Jjps  and  Aa  at  fracture,  in  addition  to  estimated  values  for  J/c  and  the  tear¬ 
ing  modulus.  Tb  develop  more  useful  expressions,  it  proves  convenient  to  define  a  common  ref¬ 
erence  position  on  the  stationary  crack  curves.  We  sdected  the  position  whmre  the  gradimit  of 
the  stationary  curve  equals  unify,  1. 

The  growth  portion  of  each  constraint  model  curve,  AJo>  is  normalized  to  define  at  the 

reference  position  using 


AX 


0  “  Atfo  +  AJpg  •  (1  —  ffljg) 


(16) 


where  m/c  denotes  the  gradient  of  the  stationary  curve  at  Jje  for  the  specific  configuration.  Fig- 
ures  14  and  15  summarize  results  for  all  the  SE(B)  analyses  (a/W=0.1, 0.5;  ns5, 10)  for  each 
of  the  Jjc  and  tearing  modulus  values  using  both  definitions  for  computation  of  the  growth 
crack  curves  (instantaneous  and  cumulative  areas).  Figure  16  provides  similar  results  for 
much  larger  selected  values  of  tearing  modulus.  Hie  normalization  using  Eq  15  accomplishes 
the  intended  purpose  of  Tninitniging  the  dependence  on  Jjc  far  growth  curves  having  a  common 
tearing  modulus.  Our  analyses  also  indicate  a  weak  dependence  of  AJq  on  ai/ao  over  the  range 
2.5:sai/ao:S3.1. 

Equations  14and  15  are  combined  to  yield  the  final  expression  for  to  correct  ex¬ 

perimentally  measured  toughness  values.  Deleting  the  notation  of  configurations  (1),  (2)  and 
the  notion  of  a  specific  point  A  from  Eq  14,  a  more  general  expression  is  given  by 

Jo  *  JS'  +  AJo  . 

(Figure  12  (a)  provides  the  graphical  interpretation  of  this  expression).  Equation  15  is  solved 
for  A  Jg  snd  the  result  substituted  into  the  above  relationship  to  define 
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^flow  ^flow 

where  normalization  of  all  J-quantities  by  the  initial  remaining  ligament  {b  =  W-  oq)  and  the 
flow  stress  is  indicated. 

Althou^  not  explicitly  used  here,  we  note  the  strong  linear  dependence  of  aJq  ^fh  AJfb 
whidi  may  lead  to  development  of  useful  expressions  of  the  form 

At/g  ■*  T^AJjpg 

where  tabular  or  curve  fit  values  for  Jk  can  be  derived  firom  Figs.  14  and  15. 

The  following  procedure  summarizes  the  application  of  Eq  17  to  correct  experimental  mea¬ 
sures  of  firacture  toughness: 

Quantitiet  needed: 

■  measured  or  inferred  values  for  n,  J/c,  Tj 

■  measured  J  at  cleavage  firacture  {Jfb) 
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FIG.  15— Groivth  portion  of  constraint  model  normalised  to  reference  point  on  stationary 
crack  solution  for  SE(B)  (alW^.S;  ns5, 10). 
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FIG.  16 — Growth  portion  of  constraint  scaling  model  normalized  to  reference  point  on 
stationary  crack  solution  for  SE(B)  (atW^O.l;  ns5, 10,  high  tearing  modulus). 
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Step*  in  procedure: 

1.  compute  A  Jre  =  Jjtb  -  Jje 

2.  look  up  AJq  using  A  in  the  graphs  provided  in  Figs.  14  and  15 

3.  look  up  m  for  Jpa  =*  «f/eusing  graphs  in  Fig.  13 

4.  evaluate  Jq  vising  the  constraint  correction  curves  for  the  stationary  crack  (  is 
Jo  evaluate  at  Jje',  Fig.  1  provides  stationary  cra^  solutions.) 

5.  evaluate  Jo  hy  substituting  thelribove  quantities  into  Eq  17 

The  constraintmodel  including  crack  growth  expressed  by  Eq  17  and  associated  figures  retains 
the  essential  features  of  original  Dodda-Anderson  modd.  This  engmeertng  modd  captures  the 
k«y  features  of  the  phenomenon  in  a  strai^tforward  approadi.  Application  of  the  new  model 
to  correct  experimental  values  of  fracture  toughness  for  constraint  involves  only  two  new  steps 
compared  to  the  model  without  crack  growth. 

Application  of  Constraint  Model  iwclniHwg  Crack  Growth 

Experimental  data  consisting  of  Jjt-Aa  values  [13,13]  that  terminate  in  cleavage  fracture 
are  presented  Figures  17  and  18  for  two  materials,  A533B  and  A515.  The  A515  data  set  [13] 
was  produced  by  testing  SE(B)  specimens  (a/W  s  0.1)  while  the  A533B  data  set  [13]  was  pro¬ 
duced  by  testing  C(T)  specimens  with  oAF  >  0.5.  The  C(T)  specimens  are  analyzed  here  using 
computational  results  for  the  deep-notdi  SE(B)  specimen.  Table  1  summarizes  the  material 
properties  and  testing  conditions.  Estimates  for  Jjc  and  for  the  constant  value  of  tearing 
modulus  indicated  in  the  table  are  found  by  fitting  the  experimental  data  (note  the  maTitniim 
amount  of  ductile  crack  extension  is  6%  of  6o  for  the  A533B  tests  and  25%  of  oo  for  the  A515 
tests).  The  rightmost  colunm  in  Table  1  shows  the  ratio  of  the  tearing  modulus  to  the  the  slope 
of  the  blunting  line  expressed  as  a  tearing  modulus.  Specific  finite  element  analyses  were  per¬ 
formed  for  these  configurations  to  accoxxunodate  the  actual  tearing  modulus  far  construction 
of  the  constraint  model. 


Sources  of  Test 
Results 

ASTM 

Steel 

Tbmp. 

[*C] 

n 

[iSpa] 

[MPa] 

[kJ/m^ 

Tj 

'^JBL 

McCabe  [13] 

A533B-G1 

-15  &0 

10 

503 

627 

75 

287 

0.35 

Kirk,etal[i3] 

A515-G70 

20 

5 

296 

545 

88 

420 

0.30 

Table  1-  Valxtea  ofJjc  and  Tj  obtained  by  fitting  experimental  toughness  values. 

Figures  17  (a)  and  18  (a)  compare  measured  toughness  values  expressed  as  Jg-ba  curves 
with  constraint  corrected  (  Jq)  values  obtained  for  each  data  point.  Filled  squares  denote  mea¬ 
sured  toughness  values  (uncorrected).  Filled  triangles  denote  Jc-values  corrected  only  for 
large-scale  yielding  effects  and  circles  denote  Je-values  corrected  for  both  large-scale  yielding 
and  cradv  growth  effects.  Jo-values  induding  growth  corrections  are  provided  for  the  model 
derived  from  the  instantaneous  stress  fidds  ahead  of  the  growing  crack  and  for  the  model 
derived  from  cumulative  areas  within  the  envdope  of  principal  stress  contours.  Jo-values 
based  on  the  cumulative  areas  are  always  larger —  they  indude  area  contributions  over  which 
the  prindpal  stress  equaled  or  exceeded  the  specified  value  during  the  history  of  growth. 

The  correction  for  large-scale  yielding  reduces  the  largest  measured  Je-values  by  approxi- 
matdy  a  factor  ofthree  for  the  deep-notch  data  set,  see  Fig.  17  (a),  and  a  factor  approaching 
five  for  the  shallow-notdi  specimens,  see  Fig.  18  (a).  Rdative  to  the  constraint  model  for  large- 
scale  yielding,  the  new  model  whidi  indudes  both  growth  and  large-scale  yielding  increases 
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A533B  (-15®,0®C),  McCabe  [18] 
C(T),  alW  =  0.52,  W  =  B  =  51  mm 

=  75kJ/m2  ^  =  400 


FIG.  11— Experimental  Jjt-6a  curves  and  three-parameter  WeibuU  distribution  for 
A533B  with  corrections  for  effects  of  crack  growth  and  large-ecale  yielding. 


the  effective  fracture  toughness  (Jq)  approximately  25-30%  for  the  deep-notdi  (A533B) 
specimens  with  almost  no  increase  due  to  growth  predicted  for  the  shallow-notdi  (A515)  speci¬ 
mens.  Ihe  A515  data  set  has  a  large  tearing  modulus  and  a  high  degree  of  strain  hardening 
whidi  combine  to  supprass  any  increased  constraint  from  crack  growth.  In  this  data  set, 
Woq/Jic  s  170  and  2.2;  the  JjpB/Job^creases  to  ">5  at  the  last  data  point.  For  compari¬ 

son,  Fig.  9  (a)  shows  the  response  of  an  SE(B)  specimen  with  a/W  =  0.1  and  n  =  5  but  with  a 
much  smaller  tearing  modulus  (TV  »  40i;s.  420 fortheA515  material)  and  smaller  deformation 
at  the  onset  of  tearing  (Wop/Jk  *  240  os.  170fortheA515  material).  In  that  case,  crack growth 
maintains  constraint  constant  at  Jfg  /Jq  ■>  2  under  increased  loading  which  implies  that  dur¬ 
ing  crack  growth  AJ/b  s  AJq;  ^  A515  data  set,  A  >  A  Jo- 

Figures  17  (b)  and  18  (b)  show  uncorrected  and  corrected  cleavage  toughness  values  in  the 
form  of  WeibuU  distributions  for  the  AS33B  and  A515  data  sets.  The  WeibuU  diagrams  show 
three  sets  of  data:  (1)  the  measured  fracture  toughness  values  (Jc),  shown  as  fiUed  squares;  (2) 
toui^mess  values  corrected  only  fn*  constraint  Qoss)  due  to  LSY,  shown  as  fiUed  triangles;  and 
(3)  tougbn«M  values  correeted  both  for  LSY  and  cradc  growth  effects,  shown  as  drdes.  The  dia¬ 
grams  iUustrate  difforences  in  the  growth  correction  derived  from  instantaneoiis  and  cumula¬ 
tive  areas  for  A533B;  for  A615  the  instantaneous  and  cumulative  corrections  are  identical  and 
only  the  cumulative  area  result  is  shown. 
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Tb  construct  the  Weibull  diagrams,  firacture  probabilities  for  the  nmked  toughness  values 
are  computed  using  Fi  =  (i-0.S)/N.  In  these  diagrams,  straight  lines  represent  the  best  fits  to 
the  distributions.  The  Weibull  slopes  of  the  dotted  straight  lines  for  the  uncorrected  distribu¬ 
tions  (2.4  and  2.6)  are  slightly  larger  than  the  theoretical  slope  of  m=2.0  for  toughness  values 
expressed  in  terms  of  J.  Similarly,  slopes  of  the  solid  straight  lines  for  the  corrected  distribu¬ 
tions  are  larger  as  well  (m  »  3.1, 3.6). 

Tire  trends  of  these  first  application  of  the  growth  correction  model  are  encouraging  with 
applicatiims  of  the  model  to  additional  dahi  sets  required  for  verification,  especially  data  sets 
with  moderate  hardening  and  lower  tearing  modulus  values  relative  to  the  blunting  line  slope. 
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SE(B),  aAV  =  0.1 ,  W  =  B  =  51  mm 

Tjj  420  ^  QQ  ^  =  700 


FIG.  18 — E:q>erimentdl  Jjt-Aa  curves  and  three-parameter  Weibull  distribution  forASlS 
with  corrections  for  effects  of  crack  growth  and  large-scale  yielding. 
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6.  CONCLUSIONS 


The  numerical  investigation  of  small  amounts  of  ductile  crack  growth  in  SE(B)  specimens 
relevant  to  cleavage  fracture  described  in  this  work  supports  the  following  conclusions: 

1 .  The  growing  crack  tip  develops  a  macroscopically  sharp  opening  profile  and  remains 
sharp  for  the  maximum  amount  of  growth  considered  in  each  combination  of  geome¬ 
try  and  material  jnnperties.  The  crack-tip  opening  angle  remains  constant  following 
the  initial  transient  period  of  growth  in  the  first  1-2  finite  elements  (Aa  a>d/c),  and 
increases  in  proportion  to  increases  in  the  tearing  modulus. 

2.  Near-tip  J-values  computed  with  a  domain  integral  formulation  approach  zero  or 
become  negative  very  near  the  growing  crack  tip.  For  domains  at  distances  firom  the 
current  tip  greater  than  0.15-0.202)o,  J  becomes  path  (domain)  independent  and 
agrees  with  J-values  computed  using  the  procedures  described  in  ASTM  E 1152-87 
(the  deformation  theory  J  derived  fiom  growth  corrected  i;  factors). 

3.  Small  amounts  of  crack  growth  can  impact  significantly  the  stress  and  deformation 
fields  ahead  of  the  extending  tip.  Compared  to  a  stationary  crack  in  an  SE(B)  loaded 
to  the  same  J,  crack  growth  can;  elevate  the  opening  mode  stresses  on  the  remaining 
ligament,  elevate  the  stress  triaxiality  ahead  of  the  tip,  and  increase  the  absolute 
size  of  principal  stress  contours  endosing  the  tip.  The  interaction  crack  growth  with 
the  degree  of  strain  hardening,  Jjc,  7j,  a/ Wand  absolute  specimen  size  is  quite  com¬ 
plex.  The  genmral  effects  of  each  parameter  are:  increased  levels  of  strain  hardening 
and  larger  tearing  modulus  values  suppress  the  effects  of  crack  growth  (drive  the 
growth  constraint  model  towards  the  stationary  cradc  model);  large  Jje  values,  small 
a/W  ratios  and  small  specimen  sizes  magnify  growth  effects  by  severely  lowering 
constraint  at  the  onset  of  growth. 

4.  The  prindpal  stress  contours  ahead  of  the  growing  crack  show  a  striking  spatial  sim¬ 
ilarity;  while  the  absolute  sizes  of  the  contoxirs  scale  nonlinearly  with  J ,  the  angular 
variation  in  shape  undergoes  negligible  changes  during  crack  growth.  These  same 
observations  about  prindpal  stress  contours  were  made  in  earlier  work  that  consid¬ 
ered  stationary  cracks  in  SE(B)  specimens. 

5.  The  JfbI^O  ratios  (FB  denotes  finite-sized  specimens)  predicted  by  the  stress-vol¬ 
ume  constraint  model  induding  crack  growth  show  a  weak  dependence  on  the  prind¬ 
pal  stress  contour  (oi/oo)  selected  for  computation.  For  SE(B)  specimens  with 
a/W=0. 1,  crack  growth  restores  a  very  weak  dependence  on  o'l/oo  relative  to  the  mudi 
greater  dependence  on  oi/ao  observed  for  a  stationary  crack  at  high  c^-values  (rela¬ 
tive  to  specimen  size).  For  SE(B)  specimens  with  a/W=0.5,  the  dependence  on  oi/oo 
develops  for  both  stationary  and  growing  cracks  when  the  global  bending  field  im¬ 
pinges  too  strongly  on  the  crack  tip  fields. 

6.  A  relatively  simple  normalization  procedure  reduces  the  dependence  of  the 
constraint  model  on  specific  values  of  Jje —  an  important  feature  wUch  increases  the 
utility  of  limited  numerical  solutions.  The  gradient  of  the  response  (curve), 

for  the  stationary  crack  model  at  J  =  Jje  provides  suffident  information 
to  capture  accurately  Jjc  effects  on  the  response  during  growth.  Dependence  on  the 
tearing  modulus  cannot  be  normalized;  however,  the  tearing  modiilus  effects  appear 
amenable  to  simple  curve  fitting  of  the  finite  element  results. 

7.  Cleavage  fracture  values  for  two  different  pressure  vessel  steels  have  been  corrected 
for  the  effects  of  ductile  tearing  and  large-scale  yielding  with  the  new  constraint 
model.  The  constraint  model  which  includes  only  large-scale  yielding  effects  appears 
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to  overcorrect  measured  Jc-valu^  in  A533B  that  have  some  tearing  prior  to  cleav¬ 
age  fracture  (Jo-values  are  too  small  with  only  the  LSY  correction).  In  contrast,  shal¬ 
low  notch  specimens  (oAV  =  0.1)  of  a  material  (A515)  with  high  strain  hardening 
(n  =  5)  and  a  large  tearing  modulus  (7^  =  400)  exhibit  essentially  no  increase  in 
Jo-values  predicted  by  the  growth  model. 
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